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ABSTRACT 


The effect on stresses in a cylindrical shell with a 
circular penetration subject to internal pressure has been 
investigated. The research is limited to thin, shallow, 
linearly elastic cylindrical shells; however, some compar- 
isons are made to thick shell experimental measurements. 
Results provide numerical predictions of peak stress con- 
centration factors around nonreinforced and reinforced 
penetrations in pressurized cylindrical shells. Analytical 
results are correlated with published formulas, as well as 
theoretical and experimental/results. An accuracy study is 
made of the finite element program for each of the con- 
figurations- considered important in pressure vessel tech- 
nology . 

A formula is developed to predict the peak stress con- 
centration factor (SCF ) for analysis and/or design in con- 
junction with the ASME Boiler and Pressure Vessel Code, 
Section VIII, Divisions 1 and 2. The formula is rationally 
derived to include all of the parameters that are required 
to define the various penetration configurations used in 
pressure vessel analysis, design, and construction. The 
accuracy of the empirical formula is determined by comparing 
to numerical, theoretical, and experimental data. In most 
cases, it is shown that the ASME Pressure Vessel Code SCF 
of 3.3 is extremely conservative. 
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CHAPTER I 
INTRODUCTION 

1.1 Statement of the Problem 

In a cylindrical shell weakened by a hole, the stress 
distribution caused by an internal pressure load applied to 
the shell will differ considerably from that in an unweak- 
ened shell. The maximum stress will be much larger if there 
is a circular hole in the shell than in the case where there 
Is no penetration. This conjecture is suggested immediately 
by the case of a flat plate weakened by a hole with the 
plate stretched per unit length in one direction and with 
one-half of this stretch per unit length in the other direc- 
tion. The maximum stress is 2.5 times the maximum stress 
in the solid plate. This factor (2.5) is known as the 
stress concentration factor (SCF). There is no reason to 
expect that the SCF is 2.5 for the shell. It depends on 
the geometry of the shell and the penetration: the curva- 

ture parameter of the shell, p a^ (p being the radius of 

RT a 

the hole, which is a circle in the projected shell surface, 

R is the radius of the middle surface of the cylinder, and 
T, the shell wall thickness); the ratio of the diameter of 
the ho lessor pipe to the radius of the shell; and the ratio 
of the thickness of the pipe to that of the shell. 
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most important feature of the stress state in the shell 
near the hole is that bending stresses occur, whereas in 
the unweakened shell only membrane stresses (— jp an(i 2 ^ ) 
are present. 

The loaded hole boundary condition in a pressurized 

cylindrical shell with a membrane or diaphragm over the 

hole to contain pressure only has been investigated by 

many authors . ^ lj ^ Another type of loaded hole 

boundary is the perpendicular intersection of two cylinders 

- shell and pipe. There are a few isolated numerical 

solutions and some theoretical invest iga- 

(1 2 11 12) 

tions. 3 3 3 There are many experimental results 

for both thin and thick shells containing nonreinforced 
penetrations (pipe only) through 21) an ^ vep y f ew 

for reinforced penetrations - pipe and pad ( a 15 » 17 ,22,24 ) 

(23 24) 

and pipe and pads. 5 

( 25 ) 

Also, the ASME Pressure Vessel Code' requires in a 
stress or fatigue analysis the stress concentration factor 
(SCF) to be not less than 3.3 for a " well designed penetra- 
tion " in a cylindrical shell unless positive evidence is 
available to the contrary. This evidence usually means a 
separate analysis to confirm the peak SCF, This factor is 
primarily needed to obtain peak stresses to perform a 
fatigue analysis to predict the remaining life at a pene- 
tration in the shell or to assure that the peak stress 
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around a penetration does not exceed allowable stresses. 

* 

The need for a more refined and/or a more clearly defined 
stress concentration factor became apparent in validating 

( 38 ) 

the useful life of 91 penetrations in a cylindrical shell^ 
located in a work/residential area of NASA - Langley Research 
Center. These penetrations in the shell range in size from 
1-inch to 60-inches in diameter. There are a few formulas 
in the published literature^ * * ’ for a membrane over 

the hole or a very thin penetration (pipe), but none are 
applicable for reinforced penetrations in pressurized 
cylindrical shells. 

The author and others were unsuccessful in obtaining 
any computer answers to an analytical finite element approach 
to the actual intersection curve of the shell, pipe, and 
pads boundaries. The compatibility equations for this 
actual curve, rather than a projected circular curve (un- 
covered during this study), were not acceptable to the 
computer program. Also, different coordinate systems for 
the shell and pipe input descriptions, solution vectors, 
and output notations were unsuccessful. Therefore, it was 
decided to abandon the analytical work and to magnetic 
particle examine and/or rework these penetrations rather 
than perform the analyses to obtain the refined stress 
concentration factors. This type of verification (field 
work in lieu of analysis) is not practical in all cases 




since this shell which contains 91 penetrations is one of 
1 3 600 pressure vessels (6000 pressure components) for which 
the structural integrity must be verified or validated in 
a five-year program at NASA - Langley Research Center. 

Thus, there is a need, for a positive and clear definition 
of a well designed penetration to allow use of the 3*3 SCF 
or to obtain the appropriate SCF. A "proven" formula to 
approximate the peak SCF and/or a finite element program 
to obtain a refined SCF would be invaluable in validating 
shells, pipes and/or pad(s) configurations. 

1.2 Object and Scope 

The objective of this study is to determine the effect 
on stresses in a cylindrical shell with a circular penetra- 
tion subject to internal pressure. The research is limited 
to thin, shallow, linearly elastic cylindrical’ shells; how- 
ever, some comparisons are made to thick shell experimental 
measurements. Results from this study provide numerical 
predictions of the peak stresses around nonreinforced and 
reinforced penetrations in cylindrical shells. Analytical 
results are correlated with published formulas, theoretical 
and experimental results. 

The present research also Investigates the convergence 
and accuracy “of different finite elements and mesh sizes. 
Finally, an approximate formula is developed to predict the 
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peak stress concentration factor for analysis and/or design 
in conjunction with the ASME Boiler and Pressure Vessel 
Code, Section VIII, Division 1 and 2.^ 2 ^ The formula is 
rationally derived to include all of the parameters that 
are required to define the penetrations used in pressure 
vessel technology. The accuracy of the empirical formula 
is determined by comparing to numerical, theoretical, and 
experimental data. Since limited data is available for 
reinforced penetrations, many different configurations are 
pursued to supplement the published data to provide the 
restrictions to the formula. These conf igurations are 
modeled utilizing finite elements where compatibility be- 
tween the cylindrical shell and the pipe/pad(s) are intro- 
duced through enforced constraint equations. The config- 
urations investigated are for force and/or pipe around 
both nonreinforced and reinforced circular penetrations in 
cylindrical shells subject to internal pressure. For the 
force case, the penetration is considered to be covered by 
a diaphragm or membrane that allows the hole edge to deflect 
and rotate. It also transmits the pressure force to the 
shell in the form of a uniform transverse shear stress at 
the hole edge. An automation computer program which 
punches cards for these configurations for input to a 
finite element computer program (NASA STRUCTURAL ANALYSIS 
PROGRAM - NASTRAN^ 28 ^) is utilized. 
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CHAPTER XI 
ANALYSIS 

2.1 Matrix Analysis by Finite Element Methods 

In the finite element method* it is necessary to obtain 
a characterization of the stiffness properties of each 
element In the structure and to relate end nodal displace- 
ments to the corresponding forces. This is expressed in 
the following form: 

[KHu} = {F> (1) 

where: [K] is the stiffness matrix 

{u} is the displacement vector 
{ P } is the force vector 

The process for generating a computer program of any 
structure that is composed of many finite elements is to 
first pick a set of local coordinates convenient for a 
typical element. The generalized element displacements 
are {£} and forces are {F}. The displacements {8} and 
stiffness matrix [K] are partitioned corresponding to ends 
I and j : 
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In the process of connecting elements, it Is found that 
one element’s local coordinates are not the same as those 
for another element. Therefore, a set of system coordinates 
is chosen that is convenient for a system of elements and 
the local coordinate points are numbered (points I, 2, 3, 
...). 

A systematic numbering process for the node points 
and members is chosen. The stiffness for each element 

Is calculated in local coordinates where i and j refer 
to the end points of each element. 

If the system coordinates or displacements are called 
(u), the transformation from an element’s coordinates to a 
system’s coordinates is accomplished by a transformation 
matrix, [a]. That is: 

{ 8 ). = [ a J { u } (3) 

The stiffness of the element is transformed to system 
coordinates by use of Equations (3), (2), and (1). 

[f] = [ a] T [K] [a] (4) 

Consider several elements that are connected. The 

next step is to generate the master stiffness matrix [K^]^ 

( 29 ) 

by summing all member stiffnesses in system coordinates. 

The compatibility equations for required coordinate 
points are introduced through multipoint constraint (MPC) 
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equations of the form: 

v A, u. = 0 (5) 

J o 

where : A is the coefficient 

u is the point 

j is the degree of freedom 

Thus, the stiffness matrix, force, and displacement 
vectors are modified by this degree of freedom link for 
each MPC Equation (5) at each grid point. (30) 

After the compatibility equations are satisfied, the 
boundary conditions (displacements of the structure) are 
enforced through single-point constraints (SPC). Finally, 
the system applied external forces { F } are identified 
and the equation: 

[KJr <u>R = <F)r (6) 

is solved. 

2.2 Formulation of the Problem 

Consider a flat rectangular plate of thickness T- 
containing a circular hole of radius P a with its edges 
parallel to axes Y ! , Z ! of the circular hole. A* material 
point or finite element grid point within the plate may be 
located through cylindrical coordinates (p, $, X) defined 
through 
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Y* .- p sin Z T = p cos . . (7) 

Suppose that the plane Y f Z* is now rolled into a circular 
cylindrical shell in such a way that the Z T -axis becomes a 
generator of the cylinder and Y f a circular arc denoted 
by Y (see Figure 2.2.1). If R is the radius of the 
middle surface of the shell and 0 the angle in any normal 
cross section of. the cylinder measured from the Y=0 plane 
in the positive direction of the Y-axis, then the following 
relationships are obtained to define the finite elements 
and grid points (see Figure 2.2.2): 

Y* = Y = R0 
R = R 

9 = Arc sin ( P - ) (8) 

Z 1 = Z = p cos <fc 


Consider that this same conf igurat ion - main (lower) 
shell (shown in Figure 2.2.1) with mid-surface radius 
R - is intersected by a branch (upper) shell with mid-- 
surface radius p 0 and thickness t, (p 0 - P a + t/2), see 
Figure 2.2.3* The axis of the branch shell is normal 
to the axis of the main shell. Both shells are considered 
to be infinite in length and capped at their ends. 

Finally, consider that these two circular cylindrical 
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Figure 2.2.1.- Coordinate system for a cylindrical shell with hole. 



p COS 0 



Figure 2.2.2.- Relationships used to define finite elements. 
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X 



Figure 2.2.3.- Intersecting circular cylinders (pipe and shell). 
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shells with mid-surface radii p 0 and R are complicated 
by the addition of one or two reinforcing pads with 
outside radius and thickness tp (see Figure 2.2.4). 

The location of the finite elements through grid points 
for a projected circular hole in a cylindrical shell 
with a pipe and pad(s) is governed by Equations (8). 

2.3 Introduction of Compatibility Equations 

The solution of these problems requires the matching 

of certain physical quantities (compatibility equations) 

along the intersection curve of the two shells and pad(s). 

When expressed in the cylindrical coordinate system, the 

intersection curve is of a very complex nature. Due to 

the difficulty of solving boundary value problems in which 

the boundaries are not situated on constant coordinate 

curves, the intersection curve can be approximated by ,p , 

p , and p equal to some constants, whenever p /R is 
K o p 0 

small. Figure 2.3.1, indicates the error involved in this 

approximation when the lower shell surface is developed 
( 11 ) 

onto a plane. It can be seen that the actual inter- 

o 

section curve does not differ appreciably from a circle, 
providing 

P, 

R 


eagb is 



< 


1/2 


(9) 
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Therefore, the actual hole boundary is assumed to be circu- 
lar, identical to that of a projected view of the penetra- 
tion. 

The compatibility equations for a pipe and shell con- 
figuration (Figure 2.2.3) are introduced through MPC- 
Equations (5). The six equations for each grid point for 
both the shell and the pipe at the pipe and shell juncture 

at p are as follows (A is deflection and 6 is rotation) 
o 


A pipe 

- 

. shell 

a r ' 

= 0 

A pipe 

a e 

- 

* shell 
A e 

T .shell _ 

2 9 Z 0 

&P lpe 

z 

- 

.shell 

A z 

T / • .shell n _ 

2 ( " e e . >- 

A Pipe 
0 R 

- 

.shell 

6 r - _ 

0 

fi pipe 

V 

- 

.shell 

e e 

0 

©pipe 

. 

Q shell 

0 


where the subscript is the degree of freedom of the grid 
point. For the finite element analysis, the independent 
degrees of freedom are those for the shell, and the depen- 
dent degrees of freedom are those for the pipe. The super- 
script in' Equations (10) denotes whether the degree of free- 
dom represents the pipe or shell. 



The compatibility equations for a "pipe-shell-inner 
pad" configuration (Figure 2.2. M) are as follows: 

1, Fipe-to-shell juncture (p^): Equations (10) 

2. Inner pad (ip) to shell junctures for each grid 
point at both p 0 and p^: 


A ip 

. shell 

- a r 

‘T 

.shell 
- A e 


. shell 

- A z 

•e P 

^shell „ 

" n 

ip 

e e 

■ .shell 

“ e e 

e^ p 

.shell _ 

- Vr? 


= o 


( T-r jp ) 6 shell } „ 0 

2 Z , 


(T+^lp) 

2 


e' heU = o 


(ID 


The compatibility equations for a "pipe-shell-outer pad(op)’ 1 

configuration for each grid point at locations p Q for the 

pipe and p and p for outer pad are as follows: 
a P 


pipe. 

R 

pipe 


.shell n 

- a r = 0 


.shell /T , , \ fi shell 

A 0 “ ( 2 + t op ) ®Z 


= 0 


(12a) 


.shell , /T , x fi shell _ 
A z + ( 2 + t o P ) 0 e 


A pip e _ 
z 


0 



e 


0 


e 


pipe 

R 

pipe 

0 

pipe 


0 


0 


0 


shell 

R 

shell 

0 

shell 


= 0 


= 0 


= 0 


(12a) 


op 

R 

o p 

0 

,op 


.shell n 

Ar - 0 


A shell ,T+ op N Q shell _ n 

A e ■ ( — 2~ : > 9 z , “ 0 


shell ( T+ op . g shell = Q 
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The compatibility equations for a '’pipe-shell-outer and 
inner pads” configuration are identical to Equations (11) 
and (12). Each of these equations is provided as enforced 
constraints to every shell, pipe, and pad(s) connecting 
grid point . 

2.4 Structural Analysis Computer Program 

The finite element method is a modern, computer-orient 





ed approach to the analysis of structures^ One of its 
principal advantages is its complete generality. This 
versatility makes it possible to consider arbitrary geom- 
etries, support conditions, loadings, and variations of 
material properties within the structures. The principal 
limitation is the cost of operation. The cost is incurred 
both in the time required to prepare the input data describ- 
ing the finite idealization of the structure and its loading 
and in the computer time required to obtain the solution. 

The process for generating the -complete finite element 
computer program is described in Chapter II, Section 2.1. 

The finite element computer program utilized in this- 
study is the latest NASA STRUCTURAL ANALYSIS (NASTRAN) 
version - level 15..5-1. Structural elements are provided 
for specific representation of more common types of con- 
struction including rods, beams, shear panels, and plates. 
The range of analyses that can be solved include static, 
elastic stability, and dynamic structural problems. 

NASTRAN has been specifically designed to treat large 
problems with many degrees of freedom. Computation pro- 
cedures in NASTRAN were selected to provide the maximum 
obtainable efficiency for large problems. NASTRAN uses a 
finite element model, wherein the distributed physical 
properties of a structure are represented by the elements 
Interconnecting at the grid points. Loads are applied 
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either at grid points or on the elements for which dis- 
placements are calculated. ^ 

The system displacements {u} in Equations (1), (3) s 
(5)» and (6) locate individual grid or node points which 
schematically represent the structure. The structure is 
approximated by connecting the se grid points with the proper 
elements (rods, bars, beams, and plates) which best describe 
the individual shapes and the overall configuration to be 
analyzed. In the process of connecting the elements to 
the grids, material .and geometrical properties (areas, 
moments of inertia, modulus of elasticity, Poisson’s ratio) 
for each element can be identified. By organizing all of' 
the grids, elements, and properties In the form acceptable 
to NASTRAN, or other general purpose finite element computer 
programs, the stiffness matrix [K] in Equations (1), (2), 
(H) and (6) can be generated in the computation process. 

The loads {F} that apply to each element and/or discrete 
grid points can be identified, and the cards generated. 

Triangular and quadrilateral elements with both inplane 
and bending stiffness are used In this study. The NASTRAN 
(level 15) numerical results were calculated using Langley 
Research Center’s CDC-6000 series computers. Also, HP-9810 
programmable desk top calculator was used to assist in 

interpreting the data and to automate the empirical formula 

■ ■ ' ■ ‘ 
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developed. The approximate number of degrees of freedom 
required to model the different configurations in this 
study are as follows: shell, 1500; pipe, 300; and pad, 

300. Once the displacements {u}, Equation (6), are deter- 
mined, internal element stresses are obtained. Finally, 
inside and outside surface stresses for each element in 
the structure are computed by NASTRAN from these internal 
stresses (membrane stress + bending stress), A shell, 
pipe, and two pad configuration (modeled with 2400 degrees 
of freedom) is presented in Appendix A, 

With regard to the cost of preparing the finite ele- 
ment program input data previously described, it is quite 
likely that this will exceed the cost of the computer 
operation in most cases. In any numerical computer method, 
the characterization process of a structure can be tedious 
and time-consuming. A major part of the cost of data 
preparation is spent In eliminating errors in the extensive 
tables of numbers required to describe the idealization. 

The extent of the input process can only be minimized by 
the use of automation. In order that the finite element 
method may be used effectively as a research, analysis, or 
design tool, it is essential that automatic mesh generation 
programs be developed which will define the idealizations 
of arbitrary shell geometries. Of similar importance to 
the practical use of such programs is automatic plotting 

ORIGINAL PAGE tS 
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to present the conf iguration and results in a readily 
usable format. 

A computer program to "automat ically " punch input 
cards in the format acceptable to NASTRAN was developed. 
This program generates the input for shell, pipe, and pads: 
grid, element, load, compatibility, and boundary condition 
cards. The cards punched from this program are input to 
NASTRAN for solving the inside and. outside surface element 
stresses. Many configurations we re solved (presented in 
Chapter V) in order to obtain trend data and comparative 
results. As a spin-off of the NASTRAN program, plots can 
be obtained for pictorial or. presentation purposes and as 
a debugging tool. Sample plots are shown in Figures 2.4.1, 
2.4.2, and 2.4.3* 


sagas? 










21 



Figure 2.4.2.- Nonreinforced hole in cylindrical shells. 
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CHAPTER III 
ACCURACY STUDY 

3.1 Introduction 

In the finite element analysis of any structure, a 
first requirement is the idealization of the structure. 

For example, a shell surface (shown in Figure 3*1-1) is 
divided into a system of appropriately shaped pieces. The 
individual pieces must be standardized as simple shapes 
such as triangles, rectangles, or quadrilaterals in order 
that their stiffness properties may be defined. This 
requirement imposes a certain degree of . approximation in 
Idealizing the geometry of shells: a curved boundary 

will usually be represented as a series of straight line 
segments. In general, this boundary approximation is not 
severe, and it can be reduced to any desired error limit 
by reducing the size of the elements. 

The most important approximation is in the shell 
behavior assumption itself. If the shell Is treated as a 
two-dimensional surface rather than as a three-dimensional 
solid, this implies certain assumptions and limitations. 
For example, in the KIrchhoff theory, it is assumed that 
stresses in the direction normal to the shell surface are 
small compared to membrane stresses, and lines normal 
to the surface are assumed to remain normal and unstrained 





Figure 3.1.1.- Discretized shell. 


25 


during deformation. Approximations of this type are not a 
special feature of the finite element solution but are in- 
herent in any shell theory. Another common approximation, 
in addition to the straight line segment representation of 
boundaries, is that the elements connecting grid points 
are flat surfaces or a group of several flat surfaces. The 
great advantage of this assumption is that the membrane and 
bending stiffness properties of the individual ‘ flat plates 
are uncoupled. The coupling, which is characteristic of 
shell behavior, is developed only in the assemblage of the 
flat plates into an approximation of the curved shell 

(31 ) 

surface . J ' 

3.2 Types of Finite Elements 

The nature of the finite element approximation is such 
that the analytical results generally converge toward the . 
true solution as the finite element mesh size is refined. 
Other factors in the convergence criteria are the number of 
grid points specified for each element, type of elements, 
and primarily, the number of degrees of freedom (DOF) at 
each grid. The Kirchhoff theory takes account of five DOF 
(3 translations and 2 rotations about axes tangent to the 
shell surface). Most shell elements make use of these same 
5 DOF at each grid. 

The finite elements employed in the discretized shell 



in Figure 3.1*1 are both planar triangles and quadrilaterals 

assembled from 4 planar triangles. The forces and stresses 

on these elements are shown in Figure 3.2,1^^ The mem- . 

brane stiffness of the triangular element is represented by 

the well-known constant strain triangle and shown in Figure 

3.2.2, The components of displacement s u and v, are parallel 

to the local coordinate system (element X and Y axes). The 

bending property based on cubic displacement patterns is 

given by a fully compatible plate bending element, a Clough 

(32 ) 

bending triangle. This triangle is formed 'by sub- 

dividing it into three basic bending triangles- as shown 
in Figure 3. '2. 3. The X-axis of each sub-triangle corresponds 
with an exterior edge, so that continuity of slope and 
deflection with surrounding Clough triangles is assured. 

The added grid point in the center is like the other grid 
points in that equilibrium of forces and compatibility of • 
displacements are required at the center point. In addi- 
tion, the rotations parallel to the internal boundaries at 
their midpoints, points 5, 6 and 7, are constrained to be 
continuous across the boundaries. The equations for slopes 
in the basic triangles contain quadratic and lower order , 
terms, and since the normal slopes along Interior boundaries 
are constrained to be equal at three points (both ends and 
the middle), it follows that slope continuity is satisfied 
along the whole boundary. Displacement continuity on all 





Figure 3.2.2.- Triangular membrane element. 



Figure 3.2.3.- Clough bending triangle. 
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boundaries is automatically satisfied when the displacement 
function contains only cubic and lower order terms. Thus, 
complete continuity of slope and displacement on all 
interior and exterior boundaries is assured for the Clough 
triangle. 

The Clough triangle is superimposed with a membrane 
triangle to form triangular elements with both membrane and 
bending stiffness. Therefore, this triangular element has 
5 DOF at each corner, 2 deriving from membrane displacements 
and 3 from the bending. The quadrilateral plate was 
developed to provide improved membrane strain behavior while 
retaining the basic 5 DOF per grid system. It is formed as 
4 ' planar triangles as shown in Figure 3-2. 4 plus 3.2.5, with 
the grids modeling the shell midsurface. Each triangle has 
one-half of the bending, stiffness or one-half of the 
thickness (membrane) assigned to the quadrilateral element. 
Since four points, in general, do not lie in a plane, care 
must be taken to ensure equilibrium and compatibility. 

Rather than try to define a warped surface, an averaging 
process is used on the noneoplanar membrane triangles. The 
bending element uses two sets of overlapping basic bending 
triangles. Since coupling between membrane stiffness and 
bending stiffness is not, at present, included in NASTRAN, . 
quadrilateral elements with both membrane and bending 
properties are treated by simple superposition of their 
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4 4 



Figure 3.2.4.- Quadrilateral membrane element. 


4 4 





Figure 3.2.5.- Quadrilateral bending element. 






31 

membrane and bending stiffness matrices. The following 
NASTRAN elements are the ones used in this research: 

1* TRIA 2 - The triangular element with bot'h inplane 

and bending stiffness. 

2. QUAD 2 - A quadrilateral element similar to TRIA 2. 
3.3 Shell 

In this section, consideration is given to the accuracy 

obtained with different mesh sizes and the two types of 

elements (TRIA 2 and QUAD 2) used to obtain the peak stress 

concentration factor for a pressurized cylindrical shell 

with a circular hole. Figure 2.2.1. A quarter of the shell 

is chosen as the model: shell radius, R; length, >_ 2R; and 

the hole radius, p . A typical finite element model is 

a 

shown in Figure 3.3.1* If the element size or general 
mesh geometry is too large at discontinuities (such as a 
penetration), the structure will be "too stiff”. The 
finite element approximations to peak- stress concentration 
factors will be below the correct answer. This model 
provides a gradual transition of large to small elements as 
the hole opening is approached. Since the "radial” section 
Is the easiest to generate or modify, different mesh 
arrangements and elements will be used to compare the 
accuracy peak stress concentration factors (SCF) to 
"exact” solutions. 


PAGINAL PAGE IS 
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Figure 3.3.1.- Typical plan view of a 
with a circular hole. 
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The rectangular arrangement of quadrilateral (HQ) 

elements of the "radial” section is shown to a larger scale 

In Figure 3-3.2. The right (RDT) and left (LDT) diagonal 

arrangements of triangular elements for the "radial” section 

are shown in Figures 3.3.3 and 3«3«^ s respectively. A 

refined version (additional grid points) of left diagonal 

triangular (RLDT) elements is shown in Figure 3.3.5= The 

triangularizing of only the first element in the rectangular 

arrangement of quadrilateral (1TRQ) elements is shown in 

Figure 3-3.6. These code names for identifying the articu- 

( 31 5 

lation arrangements are similar to those used by Melosh. 

It should be noted that efforts were pursued in refining 

these configurations by grouping the same number of grid 

points very near the hole ( see ' Figure 3 -> 3. 7 ) - all to no 

avail. Also, results for refining of the RLDT, regardless 

of the local transition, were in error. This idea of too 

( 31 ) 

much local refining is known. When the elements were 

smaller than one-half of the shell thickness, the peak 
SCF oscillated about the "exact” solution. 

A summary of the computer runs for peak stress con- 
centration factors for a pressurized cylindrical shell with 
force around the hole is presented in Table 3*3.1. The 
units used in this table and throughout this study are in 
the English or American System (inches and pounds). The. 

RQ radial section model approximates the peak inside SCF 
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to within -11# (below) of the theoretical or "exact" 

solution obtained in Reference (3)- The approximations 

for the triangular arrangements (RDT, 1TRQ, LDT) are 

converging to the theoretical solution without local 

refining (RLDT). Many authors have reported the mid- 

(12 3 ) 

surface (average) SCF. 5 5 Since this is not the 
peak value, the largest SCF (inside) will be the gauge 
in the accuracy comparisons to determine the arrange- 
ment to use for all other shell/force problems. The 
LDT (which is much easier to generate) peak inside 
SCF is slightly more accurate (-1.053S difference) 
than the RLDT (+1.14# difference). Therefore, the LDT 
arrangement is the one chosen to obtain the peak SCF 
results for shell configurations with just a force 
around the hole. 


Table 3.3-1 - Peak SCF For Finite Element Models For A 
Pressurized Cylindrical Shell With Force Around 
The Hole: p *13-0, R-112.0, v=0.3. 

a. 


RADIAL 

SECTION 

DESCRIPTION 

T 

Peak SCF <B <p= 0 and p a j 

INSIDE 

OUT 

51 D E_ 

Ava. : 1 

NASTRAN 

% DIFF 

NASTRAN 

% DIFF 

NASTRAN 

% DIFF 

RQ 

1.25 

4.997 

-10. 99 

3.108 

-16.83 

4.053 

-13.31 

RDT 

1.25 

5-405 

- 3-72 

3-395 

- 9-15 

4 . 400 

- 5.88 

ITRQ 

1.25 

5.479 

- 2.41 

3-477 

- 6.96 

4.478 

- 4.21 

LDT 

1.25 

5.555 

-1.05 

3.554 

- 4.90 

4.555 

-2.57 

RLDT 

1.25 

5.678 

+ 1.14 

3.680 

- 1.53 

4.679 

+ 0.09 

"EXACT” ^ 

5.614 

"EXACTS 

f 3.73?1 

"EXACTfi® 

f 4 .'67 5 


REF. (3) REF. (3) REFERENCES 

(1,2,3) 


JAG®® 
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3.4 Shell and Pipe 


The next configuration to be considered is two inter- 
secting circular cylinders (shell and pipe) as shown in 
Figure 2.2,3. The same shell models in Section 3. 3 are 
used; therefore, a quarter of the pipe is modeled: hole 

radius, p ; pipe radius, p ; pipe thickness^ t; pipe length, 
a o 

> 2 p , and all of the shell parameters previously described, 
— a 

Typical finite element models for a pipe are arranged and 
shown as follows: rectangular quadrilateral (RQ) S Figure 

3.4.1; right diagonal triangles (RDT), Figure 3-4,2; and 
left diagonal triangles (LDT) , Figure 3=4,3. The success- 
ful implementation of the compatibility equations for ' 
the pipe/shell juncture was accomplished by using the 
following: 

1. The 0 relationship in Equations (8) and not 
0 = 2 Arc sin ( p s . ln . , - l ) or 


0 


2 Arc sin 


2R 

, P 0 5in 
2R ~ 


5 


and 


2. The same coordinate system for both the shell 
and the pipe, and 

3. The same coordinate system for input, solution, 
and output. 

A summary of the peak SCF computer solutions for shell 
pipe is presented in. Table 3=4.1, The finite element 
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approximations are converging to the theoretical 
solutions. 3J/ For T = 1.084, the convergence is as good 
as the shell problem. For the same reasons, the left 
diagonal triangular shell and the right diagonal pipe is 
the configuration used in the numerical NASTRAN results. 

For comparison purposes, computer results for two other 
shell thicknesses (T = 1.25 and 2.167) for the "simplest" 
and the chosen configurations are also shown. 

Figures 3-4.4 through 3-4.7 show the comparison of 

the NASTRAN results with experimental and analytical data. 

The difference between the chosen NASTRAN configuration 

( 14 ) 

and the experimental data from that in Cranch is only. 

8 . 38 $, see Figure 3-4.4. Another example of. trend and 
accuracy is the nozzle-to-eylinder intersection model in 
Figure 3-4.5- The analysis was performed by Prince and the 
experimental investigation by the Oak Ridge National 
Laboratory. ' Remarkable accuracy is achieved (-3-9.5 % 
difference) with a direct trend comparison, see Figures 
3.4.6 and 3-4.7- Other examples of experimental data and 
numerical/analytical findings are presented in the results. 
Chapter V. 




i. 

























SHELL THICKNESS = 1.084 

SHELL 


INSIDE SCF 
@ <j> = 0 

ON SHELL 
& P a 

RADIAL 

SECTION 

DESCRIPTION 

PIPE 

DESCRIPTION 

NASTRAN 

% 

DIFFERENCE 

LDT 

(First Row 
Removed ) 

LDT 

2.189 

-12.09 

LDT 

(First Row 
Removed ) 

RDT 

2.197 

-11.77 


RDT 

RDT 

2. 282 

LDT 

RQ 

2. 428 

LDT 

RDT 

2.447 

LDT 

T of First 
Row = T/2 

RDT 

2.462 

RLDT 

T of First 
Row = T/2 

RDT 

2.524 


.35 


- 2.49 


- 1.73 


- 1.12 


+ 1.37 


" EXACT 
Ref (1 ) 


2.490 


SHELL THICKNESS = 1.25 


L 

T of 
Row= 


RQ 

2,213 

RDT 

2.563 

— ■ — ■ »— * — 


" EXACT” „ 
REF (1,3) 


SHELL THICKNESS * 2.167 


2.793 


■ 16 . 80 


- 3.65 


2,660 


■ 18.14 


■ 9.02 


3. 


Table 3.4.1 - Peak SCF For 
Cylindrical Shell With 
v = 0.3, And 

Finite Element Models For a 
Pipe: p ~ 13.0, R - 112.0 
t = 1.3. a 
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Figure 3^4.4. - Comparison of NASTRAN results with experimental data for shell and pipe. 










Figure 3.4.6.- Comparison of NASTRAN results with analytical and experimental 
data for shell and pipe. 



X 


Circumferential stress (outside surface of cylinder) @ <£= 0 



Figure 3.4.7.- Comparison of NAST RAN results with analyst! cal and experimental 
data for shell and pipe. 
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3 ..5 Shell , Pipe, and Pad 

The preceding sections provided comparisons of the 

accuracy of a variety of finite element solutions to 

theoretical, numerical, and experimental results. Results 

for this specific case (shell, pipe, and pad) are limited; 

however, a configuration was chosen to test the MFC cards. 

This configuration is two intersecting circular cylinders 

reinforced with a pad. The same shell and pipe models in 

Sections 3.3 and 3.4 are used; consequently, a quarter of 

the pad is modeled: outside pad radius, p^; and pad 

thickness, t ; and all of the shell and pad parameters 
P 

previously described. Typical finite element models for 

a pad are similar to those of the shell "radial” section, 

see Figures 3-5.1 through 3.5-3. A summary of the peak 

SCF computer solutions for a shell/pipe/pad configuration 

is presented in Table 3*5.1. 

There have been six shell/pipe/pad configurations 

experimentally Investigated . ^ 1 ^ 22 5 ^ The NASTRAN 

( i4 ) 

result as compared to first set of data is shown in 

Figure 3.5.4. The dashed line is the author’s extension 
of the available experimental curve. The descriptions of 
the two models to show convergence is presented in Table 
3.5.2. The second comparison is shown in Figure 3.5.5 
and presented in Table 3.5-3. This trend comparison is 
outstanding. The accuracies of the two, -3-85$ and -4.4235, 
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is remarkable. The convergence of these NASTRAN models for 
the two experimental cases lends authenticity to the first 
numerical case for which there is no comparison. Therefore, 
the LDT shell, RDT pipe, and LDT pad will be used as the 
finite element model for other numerical NASTRAN results. 

The other four experimental results will be presented in 
Chapter V (for comparison to the formula developed in 
Chapter IV). 




SHELL 

RADIAL 

SECTION 

DESCRIPTION 

PIPE 

DESCRIPTION 

TOP PAD 
DESCRIPTION 

SHELL 

SCF 

(INSIDE) 

4 > = 0 & p a 

RQ 

RDT 

RQ 

1.711 

RDT 

RDT 

RQ 

1.808 

RDT • 

RDT 

RDT 

1.831 

LDT 

RDT 

LDT 

1.868 


Table 3.5.1 - Peak SCF For Finite Element Models for 

Cylindrical Shells With Pipe and Pad: P a = 13*0, T = 1.25 

R = 112.0, v = 0.3, t » 1.3, t = 1.5, & P p = 25.0 





















Inside surface stresses @0=0 



Ramsey' s formula 

<3> NASTRAN computer run 

V WRC bulletin no. 60 
ref. (14) 

Shell Pad Pipe 

R = 24.312 p = 5.25 p Q = 3.172 

r 

T = 0.625 t = 0.625 p_ = 3.032 

p a 

t = 0.280 
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REF (14) 


TABLE 3-5-2 - PEAK SCF FOR FINITE ELEMENT MODELS FOR A 
PRESSURIZED CYLINDRICAL SHELL WITH PIPE AND PAD; 
p = 3.172, t = 0.280, R = 24.312, T = 0.625, P p = 5.25, 
t = 0.625, AND v = 0.3. 



REF (15) 


TABLE 3.5.3 - PEAK SCF FOR THE FINITE ELEMENT MODELS FOR 
A CYLINDRICAL SHELL WITH PIPE AND PAD: p = 5.0785 

t = 0.593, R - 19.0, T = 2.0, P p = 9 . 482 ,°t = 2.0, 

pags IS * p 
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3.6 Shell, Pipe, and Pads 

The shell, pipe, and pads configuration is often used 

in the design and construction of pressure vessels. ■ One 

experimental investigation for this configuration was 

(23) 

located through a paper by Kitching and Perkins. The 

experimental results were for two intersecting circular . 

cylinders (pipe and shell) reinforced with two pads (inside 

and outside of the shell) as shown in Figure 3.6.1. The 

experimental result was obtained during an investigation 

(24) 

by the British Welding Research Association (BWRA). 

The same shell, pipe, arid pad models in previous sections 

are used; moreover, a quarter of the second pad is modeled 

identical to Figures 3*5.1 and 3.5-3. A summary of the 

peak SCF computer solutions for the simplest and most . ■ 

accurate type (from previous sections) finite element 

configurations, is shown in Table 3*6.1. This is the only 

shell/pipe/pads experimental or numerical result available 

and the LDT type model has a difference of 9 % above the 

exact answer. "However, the actual maximum value (SCF) 

would be slightly higher than the one measured in the test 

since it is impossible to measure immediately at a point 

of discontinuity.."^ 2 ^ The discontinuity point would be at 

the NASTRAN peak SCF location - inside shell surface at 

<j> = 0 and p . The difference would be < 9 %, which is within 
a 
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engineering requirements. Therefore, the LDT shell, RDT 
pipe, and LDT pad(s) will be .used in modeling all other 
numerical NASTRAN results. 



FIGURE 3,6.1 - DETAILS OF NOZZLE (PIPE) IN A SHELL REINFORCED 
WITH TWO PADS IN BWRA STUDIES , ( 2*0 


SHELL 
" RADIAL” 
SECTION 
DESCRIPTION 

PIPE 

DESCRIPTION 

TOP PAD. 
DESCRIPTION 

BOT. PAD ■ 
DESCRIPTION 

SCF .. 

NASTRON 

' % 

DIFF 

RQ 

RDT 

RQ 

RQ 

1.662 


LDT 

RDT 

LDT 

LDT 




"EXACT" 

| 1.80 1 


REF (24) 


TABLE 3.6.1 - PEAK SCF FOR FINITE ELEMENT MODELS FOR A 
PRESSURIZED CYLINDER WITH PIPE AND PADS: p p = 6.13, 

t = 0.46, R = 21.8125, T = 1.625, P op = 9.86, t Qp = 1,375, 
P ip = 9.86, t ip = 1.375, AND v = 0.3. 
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CHAPTER IV 
FORMULA . 


4.1 Introduction 

In practical applications one frequently encounters 
problems in which a circular cylindrical shell is submitted 
to the action of forces distributed symmetrically with 
respect to the axis of the cylinder. The stress distribu- 
tion in wind tunnels, cylindrical containers, and circular 

pipes under uniform internal pressure are examples of such 
(34) 

problems . 

For the case of circular cylindrical shells arbitrarily 
loaded, two first approximation theories are of prime 
importance — (1) Love T s first approximation theory and 
(2) its simplified version due to Donnell. The simplified 
version led to three partial differential equations in three 


displacement components. These three equations contain 
terms which higher approximation theories have shown to be 
negligible. It is therefore permissible to simplify the 


equations by omitting such terms. ^ If these terms are 
omitted (only pressure is considered) and the thickness of 
the shell is constant, these Donnell equations lead to a 
single fourth order equation in w, the radial deflection, 
for the case of axisymmetrically loaded circular cylindrical 
shells. This equation obtained by a number of 
author s ^ 9 > 3 - 3 6 ) ls 
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d w , ET vj _ 

dz^ r 2 d ° 


(13) 


with 


D = 


ET' 


12(1 - \T) 


where w is the radial deflection' 

E is Young’s Modulus of the shell 
R is the shell mid radius 
T is the shell thickness 
p is the internal pressure 
v is the Poisson’s ratio of the shell. 

For the case of unsymmetrically loaded circular cylin- 
ders 3 the linear shallow, thin shell equations may be 
readily combined into two differential equations involving 
, only the membrane stress function F and the normal displace 
ment w. The compatibility and equilibrium equations are 

A 

3z 


(14) 


4 - 
V w + 


RD 


3F _ 
“ 


JL 

D 
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where 


v 4 = 


3Z 


29 


2 2 
az^sy* 


3y 


Elimination of the function F between the two equations 
above yields an eighth order partial differential equation 
in w of the form 


V 8 w 


Et 

DR 2 


3Z 


w _ p 


(15a) 


or 

V 8 w + 6U6 4 -Mf = -2- (15b) 

3Z 

where 

.4 = ET _ 3(1 - v 2 ) 

p — 0 0 O 9 

l6R^D l6R T 


(1-3 9-11) 

Equations (15) are known as Donnell’s linear theory. 5 
4*2 Force Around Hole 

The result of a perturbation solution to Equation (15b), 

modified to include a circular hole covered by a membrane 

2 

(Figure 4.2.1), through terms of order (3P a ) is a stress 

( 2 ) 

concentration factor at the hole-shell boundary. 

SCF * + cos 2 0 + 7r($p a ) 2 (l + |* cos 2 <J>)+ . . . 

(16) 



Hole 



Figure 4.2.1.- Penetration (radius p o ) in a pressurized cylindrical 

a 

shell (mid radius R, thickness T, and Poisson T s ratio 
v) covered by a membrane. 


co 

ro 
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For the case of 4> = 0, Equation ( 1 6 ) becomes 

SCF = 2.5 + it (BP a ) 2 (17) 

where p is the hole radius, 
a 

( 4 \ 

Savin's formula for the same problem is 


SCF = 2.5 + P 2 • (18) 

Lind's equation^®) is 

SCF -1+4 Bp a (1 + |p) • (19) 

Mershon^ 2 ”^ obtained for the same problem but with a pipe 
intersecting the shell hole 

p 3n 1/2 ? 1/2 

SCF = 2.5 + (2S-) = 2.5 +(-|f-) P 0 (20) 


where p is the pipe mid radius, 
o 

Mershon* s Equation (20) is restricted to t/T ~ 0 (force 
around hole only) or t/T small - where t is pipe thick- 
ness. 

For v = 0.3, these Equations (16 through 20) reduce to 


SCF 


Van Dyke 


* 2.5 + 


2.92 p 


a 


RT 


(2la) 


qBlGWAE PAG®® 

Of POOR §PAU1» 



= 2.5 +2.3 


(21b) 


SCF 


SCF 


SCF 


Savin 


Lind 


Mershon 


RT 


1 + 2.585 


(RT) 


1/2 


2.5 + 1.414 p 


(RT) 


T72 


(21c) 

(21d) 


The first attempt at a rational analysis of a long 

pressurized cylindrical shell having a small circular hole 

( 6 ) 

and closed at its ends is due to Lur’e. Due to errors 

introduced in the boundary conditions, his results are 

2 

incorrect. The terms in Equation (16) of order (0P a ) 
were one-half the values obtained originally by Lur’e. 

/ON 

This error was also confirmed by Eringen v ; and Lekkerkerker 
Eringen, Naghdi, and Thiel^ presented a study using the 
exact solution of partial differential equations of thin, 
shallow, cylindrical shell theory. The boundary conditions 
were satisfied by use of Fourier series and the least 
square error technique through the aid of extensive numer- 
ical calculations of the force only around a circular hole 
in a cylindrical shell. A comparison of the results from 
this "exact" approach with Equation (21) is shown in. Figure 


4.2.2. 


tffiSGINAIi 

©OP®. 
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The two experimental results are those of Bought an 
and Rothwell.^ 1 ^ In both of these experiments, only the 
membrane stresses were measured. The stress concentration 
factors correspond to the ratio of the maximum membrane 
stress at the hole to the membrane stress in the shell far 
from the hole. Excellent agreement is shown for the first 
data point (Sp 0 = 0.35)- The experimental membrane 
covering the hole consisted of a very flexible, thin metal 
plug. This may have introduced slight restraints on the 
freedom of the hole edge. By measuring only the membrane 
stress and, perhaps, introducing the slight restraint, the 
peak SCF for the second experiment would be higher due to 

the larger size hole (8p a " 0,58). 

The Van Dyke Equation (21a) is the only conservative 

( 3 ) 

equation in comparison with the ’’exact' 1 results and the 

first experimental case. Van Dyke's starting series for 

small £p a given by Equation (16), is generally accurate 
a 

to a 8p of about 0.3. ^ The required engineering 

3 . 

accuracy is exceeded for 8p > 0.6; therefore, this 

A p 

equation is the one chosen to improve on in the following 


format for the force around the hole 


SCF 


Force 


2.5 + A 



( 22 ) 


o 


totrTW AL PageJ® 

flffOOB QUALfrc 
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where 2.5 is the stress concentration factor for a flat 
plate weakened by a circular hole with the plate 
stretched per unit length in one direction and with 
one-half of this stretch per unit length in the 
other direction* 


and 


with 


2 

^ is a function to increase the stress due to 
RT 

shell curvature and thickness, and hole size. 


A for values of 


(RT) 


1/2 


shown in Figure A. 2. 3. 


The values of A to go in Equation (22) can be taken 
from this curve (Figure 4.2.3) or the maximum value of 
A = 2.7 could be used. This value is used in Equation (22) 


in this report. For v - 0.3 the result is 

.2 


SCF Force = 2 ' 5 + 2 ‘ 7 


RT 


(23a) 


or for any v 


SCF Force = 2 ' 5 + 6 ' 537 


(23b). 


(3(l-v 2 ) ) 1/i< P, 


where ” 

a 


(RT) 


1/2 



2 



Either use 



Figure 4.2.3.- Values of \ in shell with force equation. 
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4.3 Pipe Around Hole 

The installation (physically , or analytically) of a 
pipe into a cylindrical shell (Figure 4.3*1) reduces the 
peak stress concentration factor at the shell-pipe juncture. 
A function or reduction factor is needed to decrease this 
stress due to a small amount of reinforcement from the 
pipe. The reduction is dependent on the amount of pipe 
reinforcement (opening and pipe thickness) not needed to 
withstand the pressure. Thus, the following form is 
assumed 

SCF_. . = SCF_ (reduction factor). (24) 

Pipe Force 

The usual pressurization of cylindrical vessels (wind 

(37 ) 

tunnels) is analogous to a suddenly applied load. The 

dynamic response of a structure due to this type load is 
a dynamic load factor of ’’one minus some term.” Therefore, 
the above reduction factor takes the form 

2 

1 - (8p Q ) (pipe to shell thickness ratio) (25) 

where p is mid surface radius of the pipe and 8 is 
o 

modified to include Poisson’s ratio of the pipe. The pipe 
to shell thickness ratio (to be of the same order as 8) 
will be taken as (t/T) 1/Z| . Equation (24) becomes 
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SCF . = SCF _ 

pipe force 


|1 - ( 6 lPo ) 2 (^ 1/4 I • ( 26 ) 


The total equation (23 and 26) developed thus far is 


SC p = (2.5 + 6.537 ( 3P a ) 2 ) I l-(e i P 0 >"(^r) J ''' " S 


2,ta/4 
(27) 


4.-4 Pipe and Reinforcing Pad Around Hole 

The addition of an inner or outer reinforcing pad to 

a pipe in a cylindrical shell is shown in Figure 4.4.1. 

The two new parameters introduced with this conf iguration 

are the outside radius (p ) and the thickness (t ) ol the 

P ^ 

inner or outer pad. The major parameters from Equation (27) 
that contribute to this configuration's influence on the 
SCF are the mid radius (R), the thickness (T) of the shell 
and the inside radius (p ), the thickness (t) of the pipe. 

d. 

Thus, a function to decrease the stress due to a reinforcing 
pad around the pipe or hole should contain p^, R> P a > t^, 
and t or T. This function should be similar to Sp 
and multiplied times a pad to shell or pipe thickness ratio 
to some power. Note that 


BP - ^11^. 


v 2 )) 1/4 


(RT) 


1/2 “ & P a } 


■ 1/2 


( 28 ) 
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Figure 4.4.1.- Pipe (mid radius p Q , thickness t, and Poisson' s ratio i^) in a 
'■-'"■O- pressurized cylindrical shell (mid radius R, thickness T and 
Poisson's ratio v) reinforced with an inner or outer pad 
(outside radius p and thickness t ). 
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The pad to shell or pipe thickness ratio (to be of the same 

1/4 

order as $) will be taken as (t /TH) 

Therefore, the reduction in peak stress concentration 
factor due to an inner or outer reinforcing pad to a pipe 
in a cylindrical shell is 


(R P*) 


1/2 


t 

^TH ; 


1/4 


(29) 


where 


TH 


t for case of a thin, thin shell where pipe 
thickness is more important than the shell 
1 T for approaching case of a thick shell 


L where shell thickness is more important 
than the pipe. 

The limits as to when each thickness should be used 
(TH=t for | > 33 and TH=T for | ^ 33) are discussed 
in Chapters V and VI under the presentation and discussion 
of results. The SCF equation becomes 


SCF . = SCF 

pad pipe 


P - - - 

(R pa-) 


1/2 C— E — ) 

^ m t r 


1/4 


TH 


(3P) 


The total equation (27 and 30) developed thus far is 


2 , t a/4 


SCF lad' = (2 - 5 + 6 - 537(e P a } )|;L " (e iP 0 ) (_ Y _) 1 


P p , 1/2 (-V ) 174 

(R pJ 1 ' 1 TH ; 


(3D 
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where 


p is outer or inner outside pad radius and p 
P 

is outer or inner pad thickness. 


4.5 Pipe and Reinforcing Pads Around Hole 

The addition of inner and outer reinforcing pads to a 
pipe in a cylindrical shell is shown in Figure 4.5.1. The 
only new parameters introduced with this configuration are 
the outside radius and thickness of the second reinforcing 
pad. Therefore, the same type of term as Equation (29) will 
apply as a reduction to Equation (31) to obtain the peak 
stress concentration factor for two pads. The general 
SCF equation becomes 


SCF = (2.5 + 6.537(6p a ) 2 )|l " (6 1 P Q ) 2 (-|-) 1/i, | 


1/2 ( > )1/4 * wtr 1/2 5 


(32) 


1/4 


where subscripts o and. i are for the outer and inner 
pads, respectively. 

One can obtain the peak SCF for any of the previously 
developed five cases (force, pipe, outer pad, inner pad, 
and outer and inner pads) by merely substituting the 
appropriate parameters in Equation (32). For example, the 
pipe problem would be solved with - t Qp = t^ = 0, 

which would result in Equation (27).’ The practical approach 
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in the construction of cylindrical shells is to generally 
reinforce around pipes; therefore, these five configurations 
are the major thrust of this research. Another combination 
of cases can be addressed, i.e. force (t = 0) and pad(s). 

The physical significance of force, no pipe, and pad(s) 
could be (1) a cap welded to a reinforcing pad to close 
off an opening and (2) glass ports for observing the 
inside of cylindrical pressure vessels (test sections). 

The accuracy of the formula is determined by comparing 
to available and applicable published numerical, theoretical, 
and experimental data and to the analytical (NASTRAN) 
results. Many different cases and examples of this formula 
were pursued in order to uncover the restrictions in the 
formula. The accuracy, results, and applications are 
presented and discussed in Chapters V and VI. 

Since the formula is compared to theoretical (Shell 

theory) and to the author T s computer (NASTRAN) results, the 

reliability of the answers from shell theory and NASTRAN in 

this study is needed. An analogy to the cylindrical shell 

problem is the example of pure bending of an infinite plate 

( l) p ) 

with a circular hole. Reissner addresses this problem 
by including the effect of transverse . shear deformations. 

It is assumed in Shell theory and NASTRAN that the stresses 
vary linearly across the thickness. Reissner obtained a 
factor -6/5 (in lieu of 1.0) in front of the transverse 
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shear deformation terms in the stress-strain equations for 
two-dimensional plate theory. The significance of this 
difference is presented by comparing to an exact (three- 
dimensional elasticity theory) analysis of the plate with 
a circular hole subject to pure bending. The effect of 
transverse shear on the peak stress concentration factor is 
negligible provided p /T is greater than about 3*0. 
Therefore, even though Reissner’s paper is concerned with 
a flat plate, the effect of transverse shear on this study 
will be negligible as long as p /T > 3.0 and R/T > 25.0. 

d 

It is permitted in NASTRAN to use any shear factor 

(43 44 ) 

in front of the transverse shear deformation terms. J 

.The factor used in this study is 1.0. Therefore, a 

slightly different SCF than the exact answer will be 

obtained. All of the theoretical and NASTRAN results 

presented in this study have a p /T ratio greater than 4.0, 

except the NASTRAN configurations where the shell thickness 

equals 8,96 ( p /T = 1 . 2 and R/T = 13). It will be shown 
a 

in Chapter V that for experimental results for p a /T <0.6 
and R/T < 15, the transverse shear effects are not 
negligible. 
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CHAPTER V 

PRESENTATION OP NUMERICAL RESULTS 
5.1 Introduction 

The finite element technique described in Section 2.1 
has been used to complete numerical analyses for the con- 
figurations previously described (Section 2 . 2 , Chapters III 
and IV): shell with force around hole, shell and pipe, and 

shell/pipe/pad ( s ) . The units used, throughout this study 
are in the English or American System (inches and pounds).. 
Except where noted in these results, Poisson's ratio is 
taken as 0.3 and Young's modulus as 29 x 10^ psi. The 
membrane stress of a cylindrical shell Is well known.- 


a 


m 



(33) 


The ratio of the largest principal stress at a point 
in question to that which would occur at that point if the 
shell were not penetrated will be called a stress concentra 
tion factor (SCP). This is defined by 

SCF = Max, stress ( 34 ) 

a 

m 

where o is defined by Equation (33). 
m 
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The numerical (NASTHAN) results were calculated using 
Langley Research Center’s CDC-6000 series computers. All 
of the NASTRAN models are composed of triangular and 
quadrilateral elements with both inplane and bending stiff- 
ness. The pressurized cylindrical shell with a circular 
penetration is modeled as follows: 

1. General shell, Figure 3.3*1 

2. .Radial, shell section, Figure 3.3.M 

3. Pipe, Figure 3 .4.2 

4 . Pad(s), Figure 3 - 5 - 3 - 

For all of the tabulated results, the peak SCF at the 
shell/hole or pipe/pad(s) juncture is given. Appendix A 
contains a representative NASTRAN run. The detail computer 
printouts of. the numerous cases fill a volume of 7 cubic 
feet; therefore, for brevity, Appendix B contains typical 
one page summaries of the NASTRAN runs presented in' this 
section and Chapter III. These summaries for shell thick- . 
nesses 0.896 and 2.215 are organized in the following 
configurations: force, pipe, pipe and pad, pipe and 

pads, force and pad, and force and pads. 


The one page summary defines parameters used to model 
thatconfiguratlon: shell mid radius (R) and thickness (T); 

hole radios (p ); pipe mid radius (p ) and thickness (t); 

V’ . a o 

outer padi’outside radius (p ) and thickness (t ); inner 

■op op- 

pad outside radius (p^ ) and thickness (t^ ); and Poisson’s 
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radio of shell (v), pipe (v ), and pad(s) The 

location of each SCF is determined by (1) the angle (4), 
see Figures 2.2.1 and 2.2.2); (2) inside surface (I/s) or 
outside surface (0/s) of shell, pipe, and pad(s); and 
(3) ratio of the grid point location (p, see Figures 2.2.1 
and 2.2.2) divided by the pipe mid radius (p Q ). 

The peak SCF used as the design criterion is the largest 
of all SCF values computed for a configuration. 

The accuracy of the NASTRAN resuits is determined by 
comparing to theoretical and experimental data. . Since 
limited data are available for. reinforced penetrations, 
many different configurations were modeled to supplement 
the published "exact" data to provide the' accuracy and 1 ' 
restrictions to the formula. 

5.2 Formula Restrictions and Conditions 

The accuracy of the formula presented in Chapter IV 
was improved by imposing certain restrictions and conditions 
based on the results in this chapter. These results were 
from experimental, theoretical, published numerical, and 
the NASTRAN analytical results. The general restrictions 
inherent in any shell theory or finite element solution 
are as follows : 

1. -t— - < (Equation (9)) 

I\ — c. 
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2 . 


3 . 


— > 10 ^ 


f (thin shell theory) 


R 


- > 10 
T ^ 


4. p /T >_ 1.2 (Transverse shear effect) 

cl “ 


There are other conditions that the author' has used to 
improve the accuracy of the formula. These conditions are 


1. 5- < 0.4: use SCF , = SCF. _ . 

T . pipe ' force , 

but retain the calculated SCF ? or ^^pad(s) 


TH = • f t for - > 33 


T for ~ < 33 
^ 1 J 


in Equation (32). 


3. If ASME ’boded" top or bottom pad 


(39) 


(approximately 


A hole < A reinf.), see Figure 5*2.1, and ■ > 10, 

> It ' 

then treat as a two pad problem with ^ p on top 

and ~ on bottom. 

(0 P Q ) 2 (i|) 1 1.9 or pipe is "ill-conditioned." 

If (3 > treat as 'boded' pad - 

condition 3 * 

If 0.9 < p o > 2 (|) 1/4 < 1.3, use 

- ( g p ) 2 (X) 1 ^ in Equation (32) in place 
® 1 O X 

of Equation (25). 



>; 

■if 



Figure 5.2.1.- Illustration of ” reinforcement " required for ASME coded top 
or bottom pad. (38) - ■ 


CO 

r\j 
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7. If t * 0 and configuration is pad(s) reinforced, 
then p p = 2 p , TH = T,'and treat' as ’boded" pad . 

8. SCP from Equations (31) and (32) is always > 1.0. 

All the formula stress concentration factors are obtained 
from a computer program which is. based on Equation (32) and 
these restrictions and conditions. 

5.3 Force, Pipe, and Pipe/Pad (s) 

Several calculations were obtained for the configuration 

of Figure 4.2.1, force around the hole. These results were 

computed to compare with theoretical ("exact’*) result s ^ 3 ^ ^ ^ 

as well as to provide check points for the formula presented 

in Chapter IV, Equation (32). Throughout this chapter, the 

formula answers will be referred to as SCF Ramsey. The 

comparisons between NASTRAN/" exact" (N-E), Ramsey/NASTRAN 

(R-N), and Ramsey/"exact " (R-E) for different shell and 

force configurations are presented in Table 5*3.1. The 

second result in this table for - = 156.5 will arbitrarily 

T 

be defined as the only membrane result. All other results 

R 

in this table are in the thin shell theory realm, - > 10. 

T 

The overall differences of the comparisons of the results 
are as follows: N-E, 17.4$ and -t 8.9#; R-N, 8.6$ and -12.8$; 

and^-E, approximately ±5$. The minus sign in this study 

V. 

always indicates .less than and the plus sign greater than. . 



f FORCE (only) PEAK SCF % <j> = 0 & p 

( ~ — = 0.1 m > 10 

^ i ON INSIDE SURFACE S T 



R 


SCF 

SCF 

SCF 

% DIFF 

p a 

T 

NASTRAN 

"EXACT" 

RAMSEY 

M-E/R-N 

R-E 

2.42 

25 

0.2 

6.497 

5.534 

5.663 . 

17.40 
-12. 34 

2.33 

7.09 

50 

0.32 

9.603 

10.4 

10.983 

- 7.66 
14.37 

5.61 

4.84 

50 

0.40 

6.199 

5.534 

5- 663 


2.33 

8.712 

90 

0.72 

5.746 

5.534 

5.663 

Bill 

2.33 

10.842 

112 

0.896 

5.406 

5.534 

5.663 

Ha 

2.33 

13.0 

112 

1.084 

5.907 

G. 027 

6.259 

IPBtl 

3.S5 

14.52 

150 

1.2 

5.404 

■ 

5.534 

5. 663 


2.33 

13-309 

112 

1.25 

5.643 

5-746 

5.916 

HBkuI 

2.96 

13.0 

112 

1.25 

5.555 

5.614 

5.760 

- 1.05 
3.69 

2.60 

12.285 

112 

1.25 

5. 374 

5.319 

5.411 

1.03 
0. 59 

1.73 

12.25 

50 

2.0 

6.033 

6.257 

6.552 

- 3.58 


13.309 

112 

2.1667 

4 . 478 

4.457 

4.471 

0.47 
- 0.16 


10.7667 

112 

2.1667 

3.884 

3.774 

3.790 

2.92 

- 2.42 


13.0 

112 

2.1667 

4.464 

4.371 

4.380 

2.13- 

- 1.88 


12.285 

M| 

2.1667 

4.331 

4.170 

4.179 

+ 3.86 

1 

- 3.51 

13.2 


2.215 

5.233 

5.353 

5.450 

- 2.24 
4.15 

n 

12.01 

112 

2.9867 

3.836 

3.639 

3.664 

5.41 

-TT.T8 

.69 

10.752 

112 

8.96 

2.694 

2.958 

2.811 

- 8.93 

5TP" 

- 4.97 


Table 5.3.1 - NASTRAN, EXACT ^ ^ and Formula Comparisons 
for Different Shell and Force Configurations. 
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the "exact" value. In Figure 5.3.1 all three curves (NASTRAN, 
"exact" and formula) are approaching the flat plate (0 = 0 ) 
value at 2.5* For 0p < 0.6 the three curves concide. For 

d 

0 P a > 0 . 6 , the NASTRAN curve is below (as expected) the "exact", 
and the formula curve is above (as desired). 

The shell and pipe configurations (Figure 4.3*1) are 
presented in Tables 5*3*2 through 5*3*6 for different types 
of comparisons for various categories. The types of com- 
parisons are as follows: NASTRAN/exact/Ramsey , Table 5.3*2; 

analytical/Ramsey, Table 5*3.4; and experimental/Ramsey, 

Tables 5- 3. 3, 5*3*5 and 5*3*'6. The categories are defined 

as the following: thin. shell, ~ >vlO; and thick shell, 

R P r . P 

rp < 10 ; thin pipe, — - _> 10 ; and thick pipe, — < 10 . 

Several NASTRAN thin shell and thin pipe configurations were 

modeled to compare to theoretical ^ 1 1 ^ ^ , numerical^', and 
(7 14 17 20 24 ) / 

experimental v 5 J J results (see Table 5.3.2(a) and 

5.3.3(a). The overall dif f erences . of the thin shell and 
pipe NASTRAN results-, from the "exact" are -14.0# to 8 . 9 .%. 

The formula differences between NASTRAN 'and 

"exact" ^1>7 ,1 1 > 1 4, 1 7 , 20, 24 ) are " approximately + 1 ^, Note- 

the configuration where v = 0.3 and 0.5. The accuracy 

P 

of the formula as compared to both NASTRAN .solutions is 

approximately 3$. 

The thin shell and thick pipe results are shown in 
Tables 5*3*2(b) and 5.3.3(b). Those in Table 5.3.2(b) are 



o N AST RAN SCF 

A Exact SCF 

Q Ramsey SCF 



Figure 5.3.1.- Effect of shell thickness and curvature, and hole size (£p a ) on peak SCF 
for shell and force (t = 0) configurations. 
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Table 5.3.2 - NASTRAN, EXACT, and Formula Comparisons 
for Different Thin Shell and Thick/Thin 
Pipe Configurations: v = 0.3 5 v = 0.3; 

except as noted . • p 


ASTRAN 


2.812 


0.896 
0.896 
1 . 081 * 


2.1667 

2 . 166 ? 

2.1667 

2.215 

8.96 


0.576 


0.717 


3.139 

2.776 

2.760 


2.812 


0 . 717 C v p = 0.5) 3-061 


1.3 


0.96 
0.6825 
1.3 , 

0.5 


0.682 
0 . 86668 
1.3 
1.1 
0.896 


2.462 


2.721 
3 . 544 - 
2.563 
4, 728 


4.211 
2.992 
2.793 
'2.80 9 
2'. 77 0 


SCF ’’EXACT" 
Ref. ( 1 , 11) 

SCF 

; RAMSEY 

2.856 

2.898 

- 

1 — 1 
OJ 

on 

2.856 

2.844 

2.856 

2.898 

- 

2.897 

- 

3.152 

2.490 

2.107 

2.856 

2.898 

3-4 20 

3.120 

2.660 

2.562 

4 . 94 

- 

4,459 

4 .471 

3-157 

3.146 • 

• 3.07 

3.158 

3.268 

3.210 

2.790 

2.811- 


% DIFF 


N-E/R-N R-E 


- 1.5 
3-0 


24.9 

- 2.8 


5.00 


3.02 

2.97 

- 1.12 


- 4.73 


1.47 

- 8.77 

- 3.68 


5.15 

■ 9.02 


13.07 

•14.04 


(a) Thin Pipe 























































TABLE 5*3*2 - (concluded) 


V 

R 

3.032 

6.0 

24.312 

21.8125 

2.475 | 
6.0 

5.0 

21 .8125 


0.625 

1.625 


0.1 

1.625 


0.28 

1.625 


.05 

1.000 


SCF 
"EXACT" 
Ref. ( ) 

SCF 

RAMSEY 

m 

3-209 

2,407 

7. 129 (7) 
3.000 (24) 

2.956 


% DIFF. 


N-E/R-N R-E 


12.99 

-7.42 


-1.47 


(a) Thin Pipe 


p a 

R 

12.285 

112 

12.285 

112 

12.01 

112 

4.125 

22.6725 

4.125 

22.6725 

5.9375| 

22.6725 


4.845 


.938 

.938 

1.438 


NASTRAN "EXACT" RAMSEY N-E/R-N 


(2. 580) ' 2.775 ' " 
1.795 


% DIFF. 


R-E 


2.345 2.647 


. 3,3 


2.918 

2.918 

3.367 


13.51 
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(b) Thick Pipe ( Pj t < 10) 























































Table 5*3.3 - Experimental and Formula Comparisons for Thin 
Shell (R/T >_ 10) and' Thin/Thick Pipe Configurations. 


p a 

R 

1 

T 

t 

Reference 

SCF 

’’EXACT 1 ’ 

SCF 

RAMSEY 

% 

DIFF 

3.06 

11.860 

0.281 

0.25 

(17) 

3.050 

2.975 

- 2.46 

0.98 

7.659 

0.153 

- 0.021 

(20) 

4.750 

4.713 

- 0.78 

6.00 

21.813 

1.625 

0.380 

(17) 

5.000 

5.242 

4.84 




(a) THIN 

PIPE 

( P°/t 1 10) 



0.97 

. 24.844- 

1.688 

0.730 

(17) 

; 3.120 

2.524 

-19*10 

0.97 

24.844 

1.688 

0.218 



3.060 

2.561 

- 16.31 

1.91 

24.844 

1.688 

1*095 



2.980 

2.589 

-13*12 

1*91 

24.844 

1.688 

0.729 



2.980 

2.622 

-12.01 

6.00 

21.813 

1.625 

1.625 



2.700 

2.407 

- 10.85 

3.72 

24.844 

1.688 

1*593 



2.910 

2.719 

- 6.56 

6.00 

21.813 

1.625 

1.000 

1 

f 

3.100 

2.956 

- 4.65 ’ 


(b) THICK PIPE (£-0<: 10) 
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slightly more accurate than the thin pipe results for both 
the formula (+ 12 %) and NASTHAN (-4%) as compared to the 
"exact" * 5 results. For the same category, it is 

shown in Table 5.3.3(h) that the formula comparisons are 

(17) 

below the experimental by 19 % • The first four config- 
urations in Table 5.3.3(b) and the last three in Table 5.3.2(b) 
have p /T < 1.1. At the end of Chapter IV, it was stated 
that transverse shear effects would be negligible provided 
p ft /T > 3.0 and R/T > 25.0. Both of these requirements 

a — 

are violated. Therefore, the low SCF from the formula is 

to be expected since the hole is too small and/or the shell . 

too thick. Therefore, the overall accuracy of the formula 

for a thin shell and thick pipe is approximately +12% by 

not including the configurations where p /T < 1.1 and 

a 

I < J-'i.r. 

The theoretical comparison of Eringen ' s 11 ^ analytical 
results to the formula for the shell and pipe configurations 
is presented In Table 5.3*4. The accuracy between the two 
of ±8% is remarkable. The cases are tabulated by per cent 
difference from lowest to highest. These are the only valid 
cases from this report. . The others violated all four of his 
theoretical requirements: 

R/T > 10, p /t > 10, $P 0 < 0.5, and p q /R £ 1/3. 

Tables 5-3-5 and 5-3*6 present numerous experi- 
mental 21, 22,135) con f igurations for thick 
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Table 5.3.4 - Formula and Analytical ^ 11 ^ Comparison for 

Shell and Pipe Configurations. 


CASE NO. 

2R/T 

po/R 

Sp o 

t/T 

rSCF-lIT 
(Ref. (11 : 

Sf F 

RAMSEY 

% DIFF 

131 

250.00 

.050 

.35900 

.1000 

3.63 

3.330 

-8.26 

104 

100.00 

. 100 

. 45400 

.1000 

4.11 

3.823 

-6.98 

102 

100 . 00 

. 100 

; 45400 

.4000 

3.31 

3.130 

-5.43 

125 

250.00 

. 025 

.18000 

.0250 

2,86 

2.709 

-5.27 

132 

250.00 

. 050 

.35900 

.0500 

3.52 

3.337 

-5.20 

124 

250.00 

.025 

. 18000 

. 0500 

2.85 

2 .708 

-5.00 

126 

250.00 

.025 

. .18000 

. 0125 

2.85 

2.710 

-4.91 

127 

250.00 

. 025 

. 18000 

. 0062 

2.85 

2.711 

-4.89 

64 

50.00 

.050 

. 16100 

. 0125 

2.80 

2.667 

-'1.75 

130 

250.00 

.050 

. 35900 

. 2000 

3.48 

3.317 

-4.68 

20 

10.00 

. 250 

.35900 

.1250 

3.42 

3.261 

-4.63 

63 

50.00 

. 050 

.16100 

.0250 

2.79 

2.665 

-4.47 

16 

10.00 

. 100 

. 14400 

. 0250 

2.75 

2.628 

-4.42 

97 

100. 00 

.050 

. 22700 

. 0500 

2.96 

2.831 

-4.37 

98 

100.00 

. 050 

.22700 

. 0250 

2.96 

2.834 

-4.25 

99 

100.00 

. 050 

. 22700 

. 0125 

2.96 

2.836 

-4.20 

21 

10.00 

. 250 

•35900 

. 0625 

3.44 

3.302 . 

-4 . 01 

32 

25. 00 

. 050 

. 11400 

. 0125 

2,69 

2.583 

-3.99 

91 

100. 00 

. 025 

. 11400 

. 0125 

2.69 

2.584 

-3.96 

92 

100.00 

. 025 

. 11400 

.0062 

2.69 

2.584 

-3 . 94 

18 

10.00 

.250 

. 35900 

.5000 

2>82 

2.711 

-3.88 

38 

25.00 

.100 

. 22700 

.0250 

2.94 

2.831 

-3.71 

103 - 

100.00 

.100 

. 45400 

.2000 

3. 94 

3.797 

-3.64 

39 

25.00 

. 250 

. 56800 

2.0000 

2.22 

2.142 

-3.51 

123 

250.00 

. 025 

. 18000 

.1000 

2.80 

2.704 

-3.42 

31 

25. 00 

.050 

. n4oo 

. 0250 


^ — 8l 

-3.33 

90 

100.00 

.025 

.11400 

.0250 

2.67 

2.583 

-3-27 

62 

50.00 

.050 

. 16100 

. 0500 

2.75 

2.662 

-3.20 

129 

250.00 

.050 

.35900 

.4000 

3.05 

2.953 

-3.17 

37 ' 

25.00 

.100 

.22700 

. 0500 

2.91 

2.824 

-2.95 

96 

100.00 

.050 

.22700 

.1000 

2.91 

2.824 

-2.95 

15 

10.00 

.100 

.14400 

. 0500 

2.70 

2.662 

-2.89 

19 

10.00 

.250 

.35900 

. 2500 

3.27 

3.183 

-2.65 ■ 

56 

50.00 

.025 

.08000 

.0125 

2.61 

2.541 

-2.63 

57 

50.00 

. 025 

. 0800 

.0062- 

2.61 

2.542 

-2.61 

119 

250.00 

. 010 

. 07200 

.0050 

2.60 

2.534 

-2.55 

120 

250.00 

. 010 

. 07200 

.0025 

2.60 

2.534 

-2.55 

105 

100.00 

.100 

.45400 

. 0500 

3.93 

3.837 

-2.38 

11 

10.00 

. 050 

. 07200 

. 0125 

2.59 

2.582 

-2. 24 

118 

250. 00 

. 010 

. 07200 

.0100 

2.59 

2.533 

-2.18 

55 

50.00 

.025 

. 08000 

. 0250 

2.59 

'2.541 

-1.91 

89 

100.00 

. 025 

.11400 

. 0500 

2.63 

2.581 

-1.86 

117 

250.00 

. 010 

. 07200 

. 0200 

2.58 

2.533 

-1.81 

30 

25.00 

.050 

.11400 

.0500 

2.62 

2.578 

-1.61- 

10 

10.00 

. 050 

.07200 

. 0250 

2.57 

2.530 

-1.54 





Table 5.3.*» (Continued) 













































Table 5.3. 5- Experimental and Formula Comparisons for Different Thick Shell 
and Fipe (R/T and p /t < 10) Configurations. 


Reference 


SCF SCF 
"EXACT" RAMSEY 


1.41 

0.91 


0.91 
0 . 56 


0.91 

1.00 


1.4l 


0.44 


1.09 


1.41 


3.67 

3-652 


3.653 

3.286 


3.656 


19.000 


3.670 


3.694 


3.578 


3.315 


0.590 

0.553 


1.011 


0.563 

2.000 


0.590 


0.568 


0.377 


1.009 


0.445 

0.281 


0.279 

0.188 


0.400 


0.17 3. 516 


0.51 ' 3.670 


0.565 


.49 


0. 

.073 

0. 

.434 

0. 

. 236 

0. 

, 112 

0. 

.227 

0. 

■ 027 


3-070 

3.360 


3-320 

3.090 


3 . 280 
2.780 


2 . 800 


2 .'28O 


4 . 110 


2.450 


2.620 

2.920 


2.915 

2.755 


2 . 921 
2.571 


2.745 


2.749 


2.235 


4 . 105 


3.104 

2.985 


2.539 


2.683 


DIFF. 

R-E 


-14.66 

- 13.10 


- 12.20 
- 10 . 84 ' 


-10.95 
- 7. .52 


- 3.35 


- 1.82 


- 1.97 


- 0.12 


2.11 

2.93 


3.63 


6.47 

































































R 

T 

t 

Reference 

SCF 

"EXACT" 

SCF 

RAMSEY 

% 

DIFF. 

4.78 

19.00 

2.0 

2.405 


3.140. 

2.444 

-22.17 

2.88 

9.50 

1.0 

wmm 


3.779 

3.216 

-14.90 

0.91 

3.652 

0.553 

Ki 

(21) 

3-360 

2.920 

-13.10 

0.92 

3.648 

0.546 

0.280 


3.360 

2.928 

-12.86 

0.91 

3-653 

0.550 

0.285 

(17) 

3.320 

2.915 

-12.20 

4.13 

22.673 

4 . 845 

0.94 


3.300 

2.919 

-11.55 

0.95 

3.653 

0.555 

0.285 

(21) 

3.320 

2 . 940 

-11.45 

0.56 

3.286 

1.011 

0.092 

(17) 

3.090 

2.755 

-10.84 

0.91 

3.656 

0.563 

0.279 

(17, 21) 

3.280 

2.921 

. -10.95 

3.72 

24.844 

4.188 

0.593 ■ 

(17) 

3.100 

2.859 

- 7.77 

r l. 4 5 

19.0 

2.0 

1.144 

(17) 

2.750 

2.547 

- 7.38 

1.0 

19.0 

2.0 

. 0.188 

(17, 135) 

2.760 

2.571 

- 6.85 

0.92 

3.650 

0.565 

0.277 

(19) 

2.900 

2.931 

1.07 

5.94 

22.673 

4.845 

1.440 

(17) 

3.300 

3.367 

2.03 

0. 91 

3. 658 

0.565 

0.227 

(17, 21) 

2.900 

2.985 

2.93 

— ■ 

0.68 

3.658 

0.565 

0.112 

(17, 21) 

3.000 

3.104 

3.47 

0.17 

3. 658 

0.565 

0.027 


2.450 

2.538 

3.59 


0.51 


3.620 


0.490 


0.369 


2.220 


2.593 


16.80 
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shells and pipes. The comparison between the formula and 
these experiments are presented to illustrate, the need to . 
impose the restrictions presented in Section 5.2. The 
differences between the formula and thick shell and pipe 
experimental results are - 22 % to + 17 %. 

The shell/pipe/pad (top and bottom) configurations 
(Figure 4.4.1) are listed in Tables 5 . 3 • 7 . through 5 - 3 - 9 - 
The same categories as shell and pipe are maintained for 
the shell/pipe/pad comparisons. The comparisons are for 
NASTRAN, formula, and six experimental results. The thin 
shell and pipe (^ and — >_ 10) cylinders with a top pad 
are presented. in Table 5.3*7. The. differences of the 
formula as compared to NASTRAN for the top pad are -15.9# 
to 11.4$. The thin shell and pipe with a bottom pad are 
listed in Table 5.3.8. The formula differences for the 
bottom pad as compared to NASTRAN and the one experimental 
result v ' are -16.1$ to +3.2$. The thick pipe (_2. < 10), 


thin/thick shell, and pad results are shown in Table 5.3.9. 

The last three configurations in Table 5*3.7 and first 
three in Table 5.3*9 have only one difference - Poisson’s 
ratio. There are 0.8$ (Table 5.3*7) and 1.0$ (Table 5-3.9) 
differences between the NASTRAN SCF for all v’s = 0.3 and 
for v = 0.3, = V p ac } ~ 0*5. Based on these six examples, 

the formula only provides a different Poisson’s ratio for 


r/ 


Table 5-3.7 - 


p a /p o 



14.52 

15.00 


13.309 

13.65 


13.0 

13.65 


12.031 
12. 50 


NASTRAN, EXPERIMENTAL, and Formula Comparison for Thin Shell 
Thin Pipe, and Top Pad Configurations. 


SCF SCF 

NASTRAN RAMSEY 


0.72 


0.896 


0.896 


1.200 


1.25 


1.25 


1.25 


0.32 


0.34 


0.576 


0.7168 


0.7168 


0.96 


0.6825 


1.3 


0.9 38 


4.25 

11.25 

8.5 


14.7-5 


0.125 

0.25 

0.25 


0.50 


19.85935 0.50 



24.999 1.50 


1 

24.999 j 1.50 


21.5 | 1.50 


1.50 


12. 25 
12.50 

50 

2.00 

0.50 

21.0 

13.309 

13-65 

112 

2.1667 

0.6825 

2-4.999 

10.7667 

11.2 

112 

2.166? 

0.86668 

18.5625 

13.0 

13.65 

112 

2.1667 

1.3 

24.999 



2. 403 ' 2.384 - 0.79 

2.889 2.976 3.01 


2.369 2.338 - 1.31 


2.4 03 


1.945 


2.154 


2.134 


1.866 


1.924 


2 . 908 


2.421 


2.396 


2.259 


1.848 


2 . 


1.795 


1.883 


2 


-15.89 


.8 


2.030 

5.51 

2.772 

- 4.68 


8.22 

2.533 

5.72 

2.479 

9.74 














































































































Table 5*3*8 - NASTRAN and Formula Comparison for Thin Shell, 
Thin Pipe, and Bottom Pad Configurations. 


p /p 
^a o 

R 

T 

t 

P ip 

h P 

SCF 

NASTRAN 

SCP 

RAMSEY 

% 

DIFF. 

^UF 

"Exact 11 
Ref. (24) 

10.842 

11.2 

112 

0.896 

0.7168 

20.25 

1.00 

1.779 

1.835 

■ 3*15 
— — ■ 

— 

13*0 

13*65 

.112 

1.25 

1.30 

25*0 

1.50 

1.824 

1.883 

3*23 

— 

13*2 

13*75 

72 • 

2.215 

1.10 

23*0 

1.50 

3.019 

2.533 

-16.10 

— 

5.9 

6.135 

21.813 

1.625 

0.469 

9.869 

2.75 


1.827 

-3.84 

1.900 


Table 5.3*9 - NASTRAN, Experimental, and Formula Comparison for 

Thick/Thin Shell, Thick Pipe, and Top Pad Configurations. 


p a /p 0 

R 




t 

'NOTES 

SCF 

SCF 

SCF 

* DIFF 



p op 

op 

NASTRAN 

"EXACT" 

RAMSEY 

N-E/R-N 

R-E 

pLUBHiia 

19 

2.0 


• 

9.482 

2.0 


2. 399 

(15) 

2.51 

2.277 

-4.42 

-5.09 

-9.28 


19 

2.0 

* — - — — 

• 0.593 

9.482 

2.0 

v 0 * 5 

2.349 

• 

2.355 

0.26 


4.782 

5.0785 

19 

2.0 

0.593 

9.482 

2-.0 

V^a' 0 - 5 

2.372 


2. 355 

r»j 

C 

\ 


2 88 








(17) 



+1.7 

3 ’.095 

9.5 

1.0 

0.43 

4 .858 

2.0 



2.07 

2.105 


2 68 

— 







(17) 

2.149 

• • 


3.095 

9.5 

1.0 

0.43 

5*437. 

l‘. 125 



. 2.22 


-3-2 . 









(22) 




2.8815 
3.0975 

9.5 

1.0 

0.432 

5.4375 

1.125 




2.295 

2.145 


-6.5 


(a) Thick Shell. 












Table 5-3-9 ( Cont inued ) 


p /p 
K a 

R 

T 

t 

P op 

S op 

SCF 

NASTRAN 

SCF 

RAMSEY 

% 

DIFF. 

12.285 

13.65 

112 

2.1667 

2.73 

24.999 

1.50 

1.914 

2.195 

14.68 

12.01 
14 . 25 

112 

2.9867 

4.48 

18 . 0 

2. 0C 

2.212 

2.246 

1.54 

HKXC 1 

112 



1-25 

2.73 • 

24 . 999 

1.50 

1-379 

1.215 

-11.89 


(b) Thin Shell 
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the shell and pipe. The overall differences between NASTRAN 
and the formula for the shell/pipe/and pad (top and bottom) 
are -16.1$ to 11.4$. The comparisons between experiments 
and the formula are -9.3#' to 1.7? differences. 

The NASTRAN/f ormula comparisons for shell, pipe, and 
pads configurations (see Figure 4.5.1) are presented in 
Table 5.3.10. All of these results are for thin shell and 
pipe except the three thick pipe configurations denoted 
by an asterisk (#) in the ’’exact’ 1 column. The "two pad" 
experimental investigation presented in Chapter III is also 
included. The configuration that is presented in Appendix A 
as a representative NASTRAN program is indicated by 
"Appendix A" in the "exact" column. Most of the SCP values 
obtained are for configurations that represent current 
pressure vessel analysis, design, and construction - i.e. 
identical inner and outer pad thicknesses and outside radii. 
The fifth through the seventh configurations in Table 5.3.10 
have different .geometries for the two pads. The accuracy, 
of the formula for these examples is approximately +6$. 

The overall accuracy Is -8.91? to +14.65? difference. 

Figures 5-3.2 and 5-3.3 are representative of the 
manner In which peak stress concentration factors (SCF) for 
a pressurized cylindrical shell are affected by (1) thick- 
ness and curvature of shell and -hole size ($P a ), and 
(2.) reinforcement configurations (pipe and pipe/pad (s )) . 


) 












































































































































__ R = 112.0 

T =2.1667 
0 p Q = 0.5632 
v - 0.3 

Inner or outer pad 

Outside radius = p 

Thickness = t = 1 
P 

Pipe 

Pipe and pad 
Pipe and pads 
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The force value In each figure was obtained from Figure 5«3*1* 

For both 3p = 0.5632 and = 0.7414,' the largest 

o o 

reduction in SCF is due to the reinforcement either offered 

from a "thick'* pipe or from a "thin” pipe and one pad. This 

was typical of all results obtained. The SCF reduction due 

to a second pad was always minor in comparison to that 

obtained from a thick pipe or pipe/pad configuration. 

Although more pronounced for the smaller 3 p q (analogons to 

smaller hole), the peak stress concentration factors for the 

p 0 t 

shell and pipe rapidly increase as — decreases. This 
is indicative of the transverse shear effects (small p Q /T). 
Extension of the pipe and pad(s) curves to the peak SCF 
ordinate at t » 0 yields the answer for the shell/force/ 
pad(s) configurations. 

5.4 Force and Pad(s) 

It was stressed in Chapter IV that the primary thrust 
of the analysis would be for the cases previously presented. 
Another type of configuration can be solved by the formula. 
This is for t = 0 and the hole reinforced by one or two 
pads. Table 5.4.1 contains comparisons of formula to NASTRAN 
results for thin, pressurized, cylindrical shells with pad(s) 
reinforcing a hole covered with a membrane (force only). 

This type of configuration is identical to Figure 4.2.1 
plus one or two reinforcing pads. The formula results for 


0 



Table 5.^.1 




















































the top pad configuration have an accuracy of approximately 
+ 9 . 7 % } whereas all of bottom pad stress concentration factors 
were 5 % to 13 % above the NASTRAN results. All of the force 
and pads formula results were below the NASTRAN answers. 

The difference range is -10.2# to - 1 , 3 % » 
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CHAPTER VI 

DISCUSSION OF RESULTS 

This study provides numerical predictions of peak stress 
concentration factors around nonreinf orced and reinforced 
penetrations . Numerical results have been correlated with 
published formulas, as well as theoretical and experimental 
results. Most of the configurations have been for thin 
shells, and all are based on linear elastic structures. 

Some thick shell NASTRAN and published results were obtained. 
An accuracy study was made of the finite element program 
for each of the conf igurat ions considered important in 
pressure vessel technology. A formula was developed to 
predict the peak stress concentration factor for analysis 
and/or design in conjunction with the ASME Pressure Vessel 
Code. The accuracy of the empirical formula is determined 
by comparing to numerical, theoretical, and experimental 

<fata. 

The peak stress concentration factors (SCF) presented 
in Chapter V have confirmed that the 3»3 limit imposed by 
the ASME Code is indeed conservative for the vast majority 
of cases. The peak SCF results due to reinforcement (pipe 
or pad(s)) around a circular penetration in cylindrical, 
pressurized shells suggest a different number. If fact, 
78 . 5 % (124 out of 158) of the shell and pipe configurations 
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were less than 3.3. Using the procedure for a well designed 
penetration presented in Ref. (25) and shown approximately 
in Figure 5.2.1, the one pad problems (32) were checked 
for ASME Code reinforcement requirements. The six "coded" 
pad configurations are as follows: Table 5*3-7 , first two 

t = 1.25 configurations (SCF * 2.134 and 1,868); Table 5.3.8, 
second and fourth configurations (SCF = 1.824 and 1.827); 
and Table 5.3.9, first T - 1.0 and T = 1.25 configurations 
(SCF = 2.07 and 1.379). All of the other 26 one pad 
configurations were not "coded" and yet had peak SCF results 
equal to or below 3,0. The two pad problems were mostly 
"coded" (by choice). The seven configurat ions in Table 5.3.10 
which are not "coded" are the first five examples and those 
with T - 1.2 and T - 4.48. For the two pad configuration 
only one of the twenty problems had a peak SCF outside of 
the range of 1.0 to 2.0. This was a "coded" pads config- 
uration (T = 2.0 and R/T = 25) with a peak SCF = 2.177. 

Based on these results, a "thick" pipe (Figure 5.3.3), 

"coded" pad, or two pad configuration would have a peak 
SCF of about 2.2 (2/3 of the 3.3 limit), whereas the shell 
and pipe configuration could achieve the "well designed 
penetration" (SCF = 3.3) definition. 

All of the peak SCF results are summarized in Table 6.1 
according to thick, thin, or membrane shell definitions. 

The NASTRAN results compared to published values, + 17 . 4 # to 





TABLE 6.1 

SCF % DIFFERENCE SUMMARY 


THIN SHELL & PIPE 


CASES 

NASTRAN 

TO 

EXACT 

RAMSEY 

TO 

NASTRAN 

FORCE 

-8.93 

TO 

17.40 

-12.84. 

TO 

8.60 

PIPE 

-14.04 

TO 

+ 8.88 

-14.42 

TO 

14 . 28 

PIPE 

& 

PAD 

- 3.85 

-16.10 

TO 

11.42 

PIPE 

& 

PADS 

9.11 

- 8.91 
TO 

14.65 

FORCE 
& TOP 
PAD 

- 

-'9.59 

TO 

9-75 

FORCE 
& BOT 
PAD 

- 

+ 5-03 
+13.11 

FORCE 

& 

PADS 

- 

-10.20 

TO 

- 1.32 


THIN SHELL & THICK PIPE MEM. SHELL & THEN PIPE 


EXACT ' NASTRAN EXACT EXACT NASTRAN EXACT 
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-15.38 -4.03 

TO TO 
12.99 -0.55 


- 3.90 -10.85 


TO 

13,24 

-11.89 

TO 

14.68 


TO 

12.88 


-7.66 14.37 


24.91 


10 

EXACT 


TO 

NASTRAN 


TO. 

EXACT 


O 




Ill 


-8.9$, are high, as expected, when the membrane category is 
approached and low for an almost thick shell. This same 
general trend is true for the shell and pipe examples. The 
thin shell and pipe NASTRAN accuracy is +8.8$ to -14.0$. 

There is not enough data to determine any trend for the pipe 
and pad(s) NASTRAN results. 

The formula results compare favorably with NASTRAN and 
the published "experimental, theoretical, and numerical” 
data. The results for individual comparisons to NASTRAN 
and "exact" are shown in Table 6.1 by categories. The 
thin shell approximate percent difference results are as 
follows: force, +4.8$; pipe, +14.4$; pipe and pad, +13.7$; 

pipe and pads, +11.7$; force and top pad, +9,7$; force 
and bottom pad, 13$; and force and pads, -10$. 

One application of this formula is to an ASME Code 

"hillside" (elliptical hole) penetration. The amplification 

2 

factor by the code for this penetration is 1+2 sin (the 
angle the axis of the pipe makes with the normal to the 
shell wall). Two experimental results were obtained.^ 

This amplification factor was multiplied times the (circular) 
peak SCP from the formula to obtain the following comparisons 

1. "Exact" SCF « 1.840 (15) vs formula SCF = 1.826 

(-0.8$ difference) for p a 4.782, R = 19.0, 

a 

T = 1.0, t = 0.593, P = 7.^38, t = 11.625 
v = 0.3, and the angle = 20°. 
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2 , ’’Exact SCF = 2 . 8 97 ^ 38 ^ vs formula SCF = 2 o 8 97 
(-9.9% difference) for p = 11.9^ 9 R = 69*185, 

d 

T = 1.995, v = 0.3, and the angle - ^5°. 



113 


CHAPTER VII 
SUMMARY 

The effect on stresses in .a cylindrical shell with a 
circular penetration subject to internal pressure has been 
investigated. A general purpose finite element (NASTRAN) 
computer program was utilized to supplement the limited 
data for reinforced penetrations . A mesh generation 
computer program was developed to '’automatically 11 punch 

input cards in the format acceptable to NASTRAN. This 

] 

program Is readily adaptable to solving a general finite 
element shell problem. The NASTRAN compatibility equations 
for a shell/pipe or a shell/pipe/pad ( s ) configuration were 
successfully implemented. This provides a quick; access to 
a detail SCF for a reinforced nozzle in a pressure vessel. 
The accuracy of the finite element model has been Investi- 
gated. Moreover, for an immediate, approximate solution 
to this complex problem, the formula may be utilized. 

The formula Is (1) more accurate than the published 
nonreinf orced penetration formulas and, (2) believed to 
be the only formula applicable for reinforced penetrations 
in pressurized cylindrical shells. The accuracy of the 
formula was determined by comparing to the numerical, 
theoretical, and experimental data. It has been shown that 
the shell and pipe configuration can easily achieve the 
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"well designed penetration" (SCF = 3-3) definition. A 
reinforced penetration (thick pipe, coded pad, or two pads) 
would have a peak SCF of about 2,2 (2/3 of the 3.3 limit). 

This formula can be used to obtain the peak stress 
concentration factor for reinforced or nonreinf orced pene- 
trations in pressurized cylindrical shells. Thus, in the 
analysis and design of a new pressure vessel, the formula 
could save (1) the time to perform a detail analysis; 

(2) the time to construct the reinforcing pad(s) which are 
not always required; and (3) the extra cost of materials, 
fabrication, and weld examination. Also, fatigue 
analyses can be performed to obtain the remaining life in 

the penetration welds of cylindrical shells. If the 
£ 

4 

remaining life is small or exhausted, then a nondestructive 
examination (NDE, i.e. magnetic particle, ultrasonic, or 
radiographic) and/or repair of these penetrations would 
be performed. 

This type of verification (analyses and partial field 
work) of structural integrity of welds is invaluable since 
it is not practical to examine and/or repair every weld 
in all pressure vessels. The economics and feasibility of 
a validation program (6000 pressure components at Langley 
Research Center within the next five years) would be 
impractical without this formula. There have been six 
tunnels (pressure vessels) for which this procedure (formula 
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and NDE) has been followed. For example , the world's first 
known "cryogenic" tunnel utilizing nitrogen as the working 
medium incorporated this, formula in obtaining the 

expected life of this facility. 
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